a meson mass and width at finite density 



o 
o 

(N 



O 



(N 



X 



M.J. Vicente Vacas and E. Oset 

Departamento de Fisica Teorica and IFIC, 
Centro Mixto Universidad de Valencia- CSIC, 
Ap. Correos 22085, E-46071 Valencia, Spain 



Abstract 



The a meson mass and width are studied at finite baryonic den- 
sity in the framework of a chiral unitary approach which successfully 
reproduces the meson meson phase shifts and generates the /o and a 
resonances in vacuum. 
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t£j 1 Introduction 

o _ 

There is mounting evidence showing the existence of a light scalar isoscalar 
| meson [TJ E3 H|- However, its nature as a genuine meson state or as a 

resonance dynamically generated in the tttt scattering is still under debate. 
A review of the current status of the discussion can be found in ref. jHj. 

In any case, it is such a broad resonance that its effects are hardly visible 
in any phase shifts or decay plots. The situation could be quite different 
at finite densities and/or temperatures where the a could become lighter 
and narrow U\ • Apart from its implications related to chiral symmetry 
restoration, this topic is of special interest due to the relevance of the / = 
J = channel in the nucleon-nucleon interaction. Any substantial change of 
this meson properties could alter the in medium nucleon-nucleon interaction 
and therefore our current understanding of the nuclear matter in heavy ion 
collision, neutron stars or even in normal nuclei if the changes occur at low 
densities. 

The density dependence of the a properties has been studied in several 
models. Bernard et al. [HJ El HHj have found a decrease of the mass as a 
function of the density using a generalized Nambu Jona-Lasinio model, until 
it converges to the mass of the pion, its chiral partner. It is noticeable, that 
in the same model the p, u and 7r mesons have almost constant masses. 

In ref. JTI], Hatsuda et al. studied the a propagator in the linear a model 
and found an enhanced and narrow spectral function near the 2n threshold 
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caused by the partial restoration of the chiral symmetry, where the m a mass 
would approach m n . The same conclusions were reached using the nonlinear 
chiral lagrangians in ref. [T2] . 

Similar results, with large enhancements in the irn amplitude around the 
2n threshold, have been found in a quite different approach by studying 
the s — wave, I = nn correlations in nuclear matter [TS1 E3 EH E3- In 
these cases the modifications of the a channel are induced by the strong 
p — wave coupling of the pions to the particle-hole (ph) and A-hole (Ah) 
nuclear excitations. It was pointed out in jT3 CHj that this attractive a 
selfenergy induced by the 7r renormalization in the nuclear medium could be 
complementary to additional s — wave renormalizations of the kind discussed 
in [TTj IT2] calling for even larger effects. 

On the experimental side, there are also several results showing strong 
medium effects in the a channel at low invariant masses in the A(ir, 2ir) 
PHI E@ ED E2 E3] and ^(7,2vr) [23 reactions. At the moment, the cleanest 
signal probably corresponds to the A(j, 2n°) reaction, which shows large 
density effects that had been predicted in both shape and size in ref. |25] . 
using a model for the irn final state interaction along the lines of the present 
work. 

Our aim in this paper is to study the a mass and width at finite densities 
in the context of the model developed in |H1 ED ESI EH1 EH ED]. These 
works, which provide a very economical and successful description of a wide 
range of hadronic phenomenology, use as input the lowest orders lagrangian 
of chiral perturbation theory [H2] arid calculate meson-meson scattering in 
a coupled channels unitary way. The a meson is not a basic ingredient on 
this approach, it is dynamically generated and appears as a pole of the im 
scattering amplitude in the second Riemann sheet. 

The nuclear medium effects on the scalar isoscalar channel were imple- 
mented in this framework in refs. pHl ESj- As in other approaches, large 
medium effects were found. Namely, the imaginary part of the tttt scattering 
amplitude showed a clear shift of strength towards low energies as the density 
increases. 

In the next section we present, for the sake of completeness, a brief de- 
scription of the model used for the hit interaction in both vacuum and dense 
medium, which is already published elsewhere. The following section de- 
scribes the method used to calculate the nn scattering amplitude in the 
complex energy plane and to search for the a pole position at finite densities. 
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Figure 1: Diagrammatic representation of the Bethe-Salpeter equation. 

2 tttt interaction 

In this work we consider only the scalar isoscalar (a) channel and follow the 
simple method of ref. [31] for 7T7T interaction in vacuum and refs. [TT)| EBj 
for the nuclear medium effects. Additional information on this and related 
approaches for different spin isospin channels can be found in refs. j^TJ |2*HI 

123 EH 

2.1 Vacuum 

The basic idea is to solve a Bethe Salpeter (BS) equation, which guaran- 
tees unitarity, matching the low energy results to the chiral perturbation 
theory (xPT) predictions. We consider two coupled channels, tctt and KK 
and neglect the rjr] channel which is not relevant at the low energies we are 
interested in. 

The BS equation is given by 



Eq. ^is a matrix integral equation which involves the two mesons one loop 
divergent integral, (see Fig.QJ, where V and T appear off shell. However, for 
this channel both functions can be factorized on shell out of the integral. The 
remaining off shell part can be absorbed by a renormalization of the coupling 
constants as it was shown in refs. jSB • Thus, the BS equation becomes 
purely algebraic and the VGT originally inside the loop integral becomes 
then the product of V, G and T, with V and T the on shell amplitudes 
independent of the integration variables, and G given by the expression 



where P is the momentum of the meson-meson system. This integral is 
regularized with a cut-off adjusted to optimize the fit to the tt — n phase 
shifts (A = 1.03 GeV). The potential V appearing in the BS equation is 
taken from the lowest order chiral Lagrangian 



T = V + VGT. 



(1) 
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where the symbol <> indicates the trace in flavour space, / is the pion decay 
constant and M are the pseudoscalar meson and mass SU(3) matrices. 
This model reproduces well phase shifts and inelasticities up to about 1.2 
GeV. The a and /o(980) resonances appear as poles of the scattering ampli- 
tude. The coupling of channels is essential to produce the /o(980) resonance, 
while the a pole is little affected by the coupling of the pions to KK p?T] . 



2.2 The nuclear medium. 

As we are mainly interested in the low energy region, which is not very 
sensitive to the kaon channels, we will only consider the nuclear medium 
effects on the pions. The main changes of the pion propagation in the nuclear 
medium come from the p-wave selfenergy, produced basically by the the 
coupling of pions to particle-hole (ph) and Delta-hole (Ah) excitations. For 
a pion of momentum q it is given by 



(m ) q 2 U(q) 

n ? = v 2 , f ' 2 (4) 



with g' = 0.7 the Landau-Migdal parameter, U(q) the Lindhard function and 
(D + F) = 1.257. The expressions for the Lindhard functions are taken from 

ref. mm. 

Thus, the in medium BS equation will include the diagrams of Fig. El 
where the solid line bubbles represent the ph and Ah excitations. 

N o / \ / \ 'T\ / \ 

(a) (b) (c) (d) 

Figure 2: Terms of the meson-meson scattering amplitude accounting for ph 
and Ah excitation. 



In fact, as it was shown in jSTj, the contact terms with the ph (Ah) 
excitations of diagrams (b)(c)(d) cancel the off shell contribution from the 
meson meson vertices in the term of Fig. Efa). Hence, we just need to 
calculate the diagrams of the free type (Fig. |TJ) and those of Fig. Efa) with 
the amplitudes factorized on shell. Therefore, at first order in baryon density, 
we are left with simple meson propagator corrections which can be readily 
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incorporated by changing the meson vacuum propagators by the in medium 
ones. 

There are some other medium corrections, like the one depicted in Fig. El 
which has been considered in several papers like fTJ ESI Ej using the linear 



sigma model. 

In our approach, we start with only pseudoscalar meson fields. As it 
was discussed previously, the sigma is generated dynamically through the 
rescattering of the pions. A possible analog to the diagram of Fig. El is given 
by the diagrams of Fig. |3J 



Figure 4: Some tadpole diagrams contributing to the meson-meson scatter- 
ing amplitude 

We evaluate the contribution of these diagrams by using the chiral La- 
grangians jHHl EDI El E2j involving the octet of baryons and the octet of 
pseudoscalar mesons. We find these terms to be proportional to fry^p — n / ~f fl n 
and thus they vanish in symmetric nuclear matter. 

Also we would need to consider diagrams of the type shown in Fig. El We 
get a negligible contribution to the the scattering in the nuclear medium 
from this diagram jHS| and for simplicity we do not include it in this work. 

The 7T7T scattering amplitude obtained using this model exhibits a strong 
shift towards low energies. In Fig. El we show the imaginary part of this 
amplitude for several densities. Quite similar results have been found using 
different models [T3J and it has been suggested that this accumulation of 
strength, close to the pion threshold, could reflect a shift of the a pole which 
would approach the mass of the pion. Other pion selfenergy contributions 




Figure 3: Tadpole a selfenergy diagram. 
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Figure 5: Higher order tadpole diagram. 



related to 2ph excitations, and thus proportional to p 2 , can be incorporated in 
the pion propagator. As we are most interested in the region of low energies 
we can take as estimation the corresponding piece of the optical potentials 
obtained from pionic atoms data, following the procedure of ref. [HI] and 
substituting in Eq. |U 

[^f) 2 U(q) -> (£±£) 2 17(9) - 4vrC y (5) 

with p the nuclear density and Cq = (0.105 + i0.096)m~ 6 . Its effects are 
small except at large densities as can be appreciated by comparing Fig. 
with Fig. 7 of ref. ^1 where this piece is included. 

3 a meson mass and width as a function of 
the baryon density 

In this section we will consider the lightest pole position of the tttt scattering 
amplitude in the S=I=0 channel in the complex energy plane, which corre- 
sponds to the a resonance. In our model, the vacuum pole position occurs 
at a mass around 500 MeV and the width is around 400 MeV. A compilation 
of values for these magnitudes obtained in several modern analyses can be 
found in ref. [IHj. See also Fig. 1 of ref. [2] and jSj. 

The analytical structure of the meson meson scattering amplitude is 
driven by the meson loops, which are given by the formula 

G™ = ij o jAiD^q\q)D^s-q Q ,-q) (6) 

where D n is the ir propagator which, in the nuclei, incorporates the p-wave 
coupling of the pions to particle-hole and A-hole excitations. We can sim- 
plify the calculation by using the Lehmann representation for the meson 
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Figure 6: Imaginary part of the tttt scattering amplitude at several densities. 



n 2 f°° x Im D n (x, q 



propagators 
After some algebraic manipulation we obtain 

where F(E) is a real function given by 

F{E) = J ±jL J E E ^ImD 7T ((E + x)/2,q)ImD w ((E-x)/2,q) (9) 
which in vacuum takes the simple form 



FUE) = = -J- - 4 (10) 

As we are interested in the positive energy region, we can concentrate our 
attention on the first term of Eq. |H] which is the only one that has a cut. We 
can rewrite that term, neglecting the e, like 

F(E) F(Real(P )) + F(E) - F(Real(P )) 



P° — E P° — E P° — E 
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The E integration of the first term it is the only piece that depends on the 
Riemann sheet and can be readily evaluated. In the second sheet the integral 
takes the value 



L, dE( P°-B = L ^ P°~E > + (12) 



Substituting this result in Eq. IHlwe get the following expression for the two 
meson propagator in the second sheet 

CUP = iF^aKP")) + jT f - ^) F(S) (13) 

The rest of the pieces entering the BS equation are analytical single valued 
functions and therefore we can proceed to look for poles in the scattering 
amplitude. 



4 Results 

The function F(E) includes the phenomenological information on the pion 
selfenergy in the nuclear medium. The results for several densities are shown 
in Fig. 

The F(E) strength is related to the imaginary part of the a propagator 
and therefore reflects the energy and density dependence of the different o 
decay channels. We can start classifying them according to their density 
dependence (See. Fig IE}. 

At low densities, the a meson can only decay into two pions, a — > n + ir. 
Therefore we have a threshold at E — 2m n . As the density grows the a — > 
7r + (Nh) decay channel becomes relevant and we have a new threshold at 
E = m n , which is clearly visible in the curve corresponding to p = po/8. 
Finally, at larger densities, mechanisms such as a — > 2(Nh) which go like p 2 
become important. They are possible even at very low energies. 

Using the two meson propagator of Eq. [H] we can solve the BS equation, 
obtain the meson-meson scattering amplitude and look for the poles in the 
complex energy plane. The results for the a pole position are shown in Fig. 
E3 for densities up to 1.5 po- Note, however, that the calculation is more 
reliable at low densities because some contributions of order p 2 or higher 
are missing. We find that both mass and width decrease as the density 
increases, reaching a mass around 250 MeV and a similar width at 1.5 times 
the nuclear density. At large densities the 2ph pieces of the pion selfenergy 
become relevant decreasing further both mass and width. These results are 
little sensitive to the kaon channel. Its omission affects less than 1 percent 
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Figure 7: F(E) at several densities. 



the mass and increases the width around 5 percent. The results of Fig. El 
can be cast in terms of an effective potential which can be approximated by 

V p = a^- + b(^-) (14) 
Po \PoJ 

with a = -358 - i 108 MeV and b = 140 + i 23.6 MeV. 

Qualitatively similar results for the mass are found in other models. See 
for instance Fig. 10 of ref. [IS!- However, we should stress that there the 
changes respond to a reduction of the pion decay constant (/) value which 
we have kept constant. The basic ingredient that drives the mass decrease 
in our calculation is the p-wave interaction of the pion with the baryons in 
the medium. 

Although the a mass drops significantly, the width stays relatively large, 
even when the mass is close or below the 2ir threshold. This width comes 
mostly from medium decay channels, namely, the decay into a single pion 
and a ph excitation or 2ph excitations. Therefore we should not expect any 
signal of narrow a mesons in the medium as expected in [Zj or in models that 
include only purely mesonic decays. Nonetheless the low mass could modify 
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0(0) 0(p) 0(p 2 ) 



Figure 8: a decay channels classified according to their density dependence 



the long and medium range nucleon nucleon interaction in nuclei or nuclear 
matter. 

5 Conclusions 

In this work, we have searched for the lightest pole of the nir scattering 
amplitude, the a meson, in the complex energy plane as a function of the 
baryon density. 

The growing consensus that the a is rather a tctt resonance of dynamical 
origin than a genuine QCD state built from qq pairs, finds its support in 
recent chiral unitary formulations of the tot interaction which we have fol- 
lowed in the present work. This approach has been extended to include the 
interaction of the pions with a nuclear medium, which allows us to trace the 
density dependence of the a pole. 

A quite general agreement has been reached about these medium mod- 
ifications of the 7T7T interaction between different groups, either working in 
the present approach or using other irir scattering formulations, as long as 
they satisfy some minimal chiral constraints. All these models are, however, 
involved numerically to the point that extrapolating analytically the results 
to the complex plane is a difficult task, particularly when some pieces of the 
7r selfenergy, like the 2p2h contributions, are taken from a comparison with 
pionic atoms data and not from a theoretical model. 

In this work, the analytical extrapolation was made possible by the use 
of the Lehmann representation for the it propagators. This allowed us to 
express the scattering matrix in terms of the pion propagators evaluated on 
the real energy axis. 

The results obtained show a dropping of the a mass as a function of the 
density, down to values close to two pion masses at normal nuclear density. 
The width decreases moderately and reaches values around 300 MeV at p = 
Po. Yet, this decrease is surprising in view of the fact that in the nuclear 
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Figure 9: a mass and half width as a function of the density. Dashed lines 
include also the 2ph pion selfenergy pieces 



medium there are more channels for the a decay, like pion-p/i or 2ph. The 
presence of these channels partly makes up for the strong reduction of the 
2tt channel due to the smaller phase space available when the a mass drops. 

Altogether, the changes in the mass and width of the a are quite large, 
compared to other typical meson or baryon properties in a dense medium. 
Therefore, we expect drastic signals in the invariant mass distributions of two 
pions (or photons) in 2% (27) production experiments if those pions (photons) 
can couple strongly to the scalar isoscalar channel. Some present data seem to 
support these claims, although further experiments and calculations should 
be performed to test these findings. 
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